Generalized Quantum Walk 
in Momentum Space 



A. Romanelli, A. Auyuanet , R. Siri, G. Abal, 
and R. Donangelo ^ 

Instituto de Fisica, Facultad de Ingenierta 

Universidad de la Repilblica 
C.C. 30, CP. 11000, Montevideo, Uruguay 



Abstract 

We consider a new model of quantum walk on a one dimensional momentum space 
that includes both discrete jumps and continuous drift. Its time evolution has two 
stages; a Markov diffusion followed by localized dynamics. As in the well known 
quantum kicked rotor, this model can be mapped into a localized one-dimensional 
Anderson model. For exceptional (rational) values of its scale parameter, the system 
exhibits resonant behavior and reduces to the usual discrete time quantum walk on 
the line. 
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1 Introduction 

The quantum random walk on the line has been studied as a natural gener- 
alization of the classical random walk in relation with quantum information 
processing. Two cases have been considered: the quantum walk with continu- 
ous time or with discrete time steps [1]. One of the most striking properties of 
the quantum random walk on the line is its ability to spread over the line lin- 
early in time, while its classical analog spreads out as the square root of time. 
Experimental schemes to implement the quantum walk have been proposed 
by a number of authors [2,3,4,5,6]. 
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Another subject that has drawn much attention in qTiantiim information pro- 
cessing is dynamical locahzation (DL) [7,8,9,10], due to its potential for process 
control in a quantum computer. DL is a quantum interference effect which sup- 
presses quantum diffusion. It is typical in one dimensional periodically driven 
systems with chaotic classical counterparts, such as the quantum kicked rotor 
model [11,12]. In recent works [13,14], quantum algorithms which simulate the 
quantum kicked rotor faster than classical algorithms were presented and we 
have proposed the use of these algorithms to describe the evolution of discrete 
time quantum random walk [15]. 

In this work we consider, in section 2, a generalization of the quantum random 
walk on a line, along a line similar to the stroboscopic discrete-time quantum 
walk presented in [16,17], is considered. In section 3, we numerically character- 
ize its dynamical properties focusing in its localized and resonant behaviors. In 
section 4, we establish a connection with the one-dimensional Anderson model 
from solid state physics and discuss similarities of the modified quantum walk 
with the kicked rotor model. The last section contains our conclusions. 



2 Modified quantum walk 

The discrete-time quantum walk on the line describes a particle which is free 
to move over a lattice of interconnected sites labeled by an index k. In the 
classical random walk, a coin flip randomly selects the direction of the motion. 
In the quantum walk, the direction of the motion is selected by introducing an 
additional degree of freedom, which we call the chirality, which can take two 
values: "left" or "right", \L) or \R), respectively. At each time step, a rotation 
(or, more generally, a unitary transformation) of the chirality takes place and 
the particle moves according to its final chirahty state. The global Hilbert 
space of the system is the tensor product He where ifg is associated to 
the motion on the line and He is the chirality Hilbert space. 

The conditional translation can be written as S'^'^ where S is the one-site 
translation operator, S\k) = \k -\- 1), and az the usual 2x2 Pauli matrix. 
Then, this operator shifts the position on the lattice left or right, according 
to the chirality component 

s'^^mk) = \L,k+i) 

S''^\R,k) ^\R,k-l). (1) 

In this work, the unitary operation on chirality is a Hadamard operation, 
H — {ax + (Tz) I -\/2. Then, the Hadamard walk on the fine evolves in a time 
step T with U = S'^'^H as 

\^{t^T)) = um))- (2) 
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The sites \k) arc usually associated with position cigcnstatcs, but they can 
equally well be interpreted as momentum eigenstates, i.e. P\k) = hk\k), lead- 
ing to a quantum walk in momentum space. In fact, this is actually the case 
in proposed implementations of the quantum walk using classical waves [5,6]. 
Interpreting the lattice sites as momentum eigenstates allows us to establish 
a connection with the kicked rotator. 

We modify the quantum walk described in eq. (2) by considering that between 
two consecutive applications of U the random walker drifts with constant 
momentum. This drift is represented by a site-dependent phase shift related 
to the kinetic energy P^/2m, so that the new evolution takes the form 

\^U{t + r)) = e-'^^^^'^'/^'Ul^it)), (3) 

where f2 = and m is the particle mass. The wave vector can be expressed 
as a spinor 




with the upper (lower) associated to the left (right) chirality. The unitary 
evolution implied by eq.(3) can be written as the map 

a^it + r) = ^ [ak+i{t) + 6,+i(t) ] e'^^'^^'^' 

hit + r) = ^ [ak-i{t) - h^,{t)] e-''-'''\ (5) 

Note that for integer values of il the usual Hadamard walk is recovered. 

The evolution of the momentum distribution, Fk(t) = |afc(t)|^ + \bk{t)f, is 
obtained from eqs.(5) as in our previous work [15,18], 

Fk{t + T)^^ [Fk+i{t) + Fk-i{t)] + Pk+i{t) - Pk-i{t), (6) 

where Pkit) = 3^ [<^fc(^)^ik(^)]) is an interference term with ^{z) indicating the 
real part of z. If these terms are neglected in eq.(6), the resulting evolution 
is markovian. In the continuum limit, taking r — 1, the familiar diffusion 
equation results 

dF ^ 16^ 

dt 2dp^' ^ ^ 

The dynamics defined by eqs. (5) presents a very rich behavior, from quan- 
tum diffusion to quantum localization. For special values of a resonant 
dynamics results, in which the variance of the distribution cr"^ = J2k k'^^kit) — 
[J2k k Fk{t)]^ , increases quadratically in time. In fact, as mentioned before, 
for integer ft the usual quantum walk, eq. (2) is obtained. As is typical of 
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Fig. 1. Variance o"^ as a function of time for several irrational values of the parameter 
ri, namely 2Tri^ = 0.1,0.2,0.3 (thin lines). The walker has been initialized at the 
origin in the state, |^'(0)) = ^(l,i)|0). The dashed line represents the classical 
linear diffusion. The approximately parabolic full line corresponds to the rational 
case Q = 1/9. 

near-resonant systems, small variations in the parameter Q may lead to very 
different dynamics, as shown in Fig. 1. 



3 Localization and Resonances 



Two clear cut different behaviors are found depending on the rational or irra- 
tional value of Q. When Q is irrational, an initial stage of quantum diffusion 
with the diffusion coefficient expected from eq. (7) is followed, after a char- 
acteristic time, by a localized dynamics in which the growth ceases due to 
quantum interference (sec Fig. 1). The exponential character of this localiza- 
tion is apparent from Fig. 2, where the distribution Fk{t) is shown after 2000 
time steps, when DL has already set in. 

In Fig. 3 the loss of localization that takes place for rational values of Q is 
illustrated by considering a set of near-resonant parameter Q = Qr + S with 
Qr rational and 5 <^ 1. For not too small values of 6 {fl approximately irra- 
tional) the distribution is exponentially localized. As 5 is gradually decreased 
the exponential localization transforms into an almost uniform distribution. 
To aid visualization we have plotted only the non-zero probabilities in the dis- 
tributions since the map (5) implies that if a single site is initially occupied, at 
any later time only odd or even sites are occupied with non-zero probabilities. 
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Fig. 2. Momentum probability distribution Fk{t) at t = 2000, as a function of the 

momentum index k for 217^1 = 0.1,0.7. The dashed hnc is the exponential fit for 
2ttO, = 0.1 from where it is possible to obtain the localization length. The initial 
condition is the same as in Fig. 1. Notice that the vertical scale is logarithmic. 
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Fig. 3. Momentum probability distribution Fk{t) at t = 2000, in logarithmic scale, 
as a function of the momentum index k iov ^ = ji+S , with 5 = 10~^, 10~^, 10~^and 
10~^. The dashed line is the exponential fit for S = 10~^. The initial condition is 
the same as in Fig. 1. 

When Q is rational, the variance increases quadratically in time, as shown 
in Figs. 1 and 4. In this case, the characteristic transition frequency between 
momentum eigenstates is commensurable with the random walk's frequency. 
This is the condition for resonant behavior. To understand such a resonant 
regime, it is possible to use the theoretical developments presented either in 
our recent work [15] where we consider the principal resonances {Q integer) of 
the quantum random walk or, alternatively, in the detailed results obtained 
for a related system, the kicked rotor, in the past [12]. If Q is a non-integer 



5 



80000 



60000 



40000 



1 ' 


1 ' 1 


/ / 


- 




1/3 / 






- 




3/5 














— T 1 







Qi rr-r^— -^-r^ 1 1 \ ' L 

200 400 600 800 



Time 

Fig. 4. Quadratic growth of the variance as a function of time for several rational 
values of the parameter Q = 1/3, 1/5,3/5, 1/7, 1/9. 



rational number = |, we are in the presence of secondary resonances. The 
variance still displays a quadratic growth but it is slower than the growth 
associated to the main resonances. This is due to the fact that while in the 
main resonances all phases add constructively, in secondary resonances there 
are (at most q) different phases that repeat themselves periodically in time. 
These phases determine the time required for coherent behavior to influence 
the dynamics. 
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Fig. 5. In the upper frame, momentum probability distribution F]^ as a function 
of the momentum k for O = ^. In the lower frame the same distribution for 
O = + 10~^. In both frames the distribution is calculated for t=2000r. 

In Fig. 5 (a), the probability distribution after 2000 steps of the secondary 
resonance = shown. A constructive interference phenomenon can be 
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clearly seen when the momenta are multiples of the denominator g = 11, 
i.e. k = —22, —11, 0, 11, 22, .... On the other hand. Fig. 5 (b) shows that 
the distribution for = ^ + 10"*^, after the same time number of steps, is 
clearly localized. Thus there is a high sensitivity of this system to small varia- 
tions of Q near a rational value. Although Q is strictly rational in both cases, 
for the second, the time at which the system becomes sensitive to its rational 
character is much greater than the time considered in the figure. 



4 Connection with Anderson localization 

As we have shown, the map (5) generates the same qualitative behavior as that 
of the quantum kicked rotor. This last system has been extensively studied 
because its evolution operator can be expressed analytically and its classi- 
cal version (the standard or Chirikov map) is a cornerstone in the theory of 
chaotic Hamiltonian systems [19]. In particular, the phenomena of DL was 
first observed in numerical experiments [11] for the quantum kicked rotor. 
In the last decade, direct observations of DL have been done using samples 
of cold atoms interacting with far-detuned light [20]. When the light field is 
switched on and off periodically, the resulting Hamiltonian is essentially that 
of the quantum kicked rotor [21,22,23]. In a particularly interesting theoretical 
development [24,25], the kicked rotor Hamiltonian has been mapped into an 
Anderson model in a one-dimensional lattice [26,27]. This connected the dy- 
namical exponential localization in momentum space to the spatially localized 
wave functions of interest in solid state physics. 

There is a strong similarity between the modified quantum walk and the kicked 
rotor. The Floquet {i.e. one-step evolution) operator of the kicked rotor in- 
volves a kick operator e^*^'^°^^, with K determining the kick strength and 6 
the angular position of the rotor, followed by a free evolution over a time step, 
g-ip2/2?i p^j. ^YiQ modified quantum walk, the kick operator is replaced by the 
conditional translation and the Hadamard rotation U — S'^^'H. 

Not surprisingly, the modified quantum walk on the line can also be mapped 
into a localized one dimensional Anderson equation. Thus, this momentum 
localization may also be thought of as a dynamical form of Anderson localiza- 
tion. In order to show this, we relate the Floquet eigenstates of the modified 
walk with the exponentially localized states of the Anderson model. 

The characteristics of the modified quantum walk imply that the Hamilto- 
nian of the system is periodic in time and Floquet theory [28,29] applies. 
Consider a Floquet eigenstate \^w) of the one-step operator e^*^'^^^ U, so 
that e^*^'^^-^^/'^^ Ul'^w) = e^*™|\I'^) and w is the associated quasienergy. Then 
1^^) = e-^"'*|$^(t)) with + r)) = |$^(t)). The time dependence of Flo- 
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quet states is trivial and wc shall suppress it from our notation. In the case of 
Floquet eigenstates, the map (5) reduces to the set of equations 



fk^k — flfc+l + bk+l 

fkbk = Ofc-i — bk-1, (8) 

where fk = y^e**-^'^^*^^"'^-'. Note the dependence on the quasienergy w implicit 
in (8). After decoupling these equations and performing the change of variable 



(^k\ -k p I 



we obtain 



hPk = Pk+1 +Pk-1: (10) 

where we have defined gk = i{fk+i - fk) and gt = i{fk-i - fk)- 

These equations strongly suggest a relation with Lloyd's model [30], a tight 
binding model of an electron on a one-dimensional disordered lattice. This 

model is a special case of the one-dimensional Anderson model [27] in the 
case of nearest neighbor interaction, which has been analytically solved [29]. 
The difference with our model is that gk and gk in eqs. (10) have a non-zero 
imaginary part. 

It is clear that the equations for both spinor components have same nature, 
so we focus on the first and show that it can be transformed into an Anderson 
equation 

nak + J^Wkm-^^ (11) 

Ijtk 

where Tk is the kinetic energy at site k and Wki is the hopping term. It is 
well established [29] that if the are pseudo-random functions of the site k 
and the hopping term Wkm decays sufficiently fast with distance \k — m\, the 
wavefunction is exponentially localized. We use eq. (10) recursively to express 
it as an equation with real coefficients of the form (11) with 

Tk = gi [gi-i + gUi] + gUA-i [{gl)' + (gi)'] , (12) 

gl and gl are the real and imaginary parts of gk- 
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The potential term Wki is 



-9i-2 

'9l+2 



9l-2 9l-3 

9i9i-i + 9i9i-i 
9i9i+i + 9i9i+i 

9l+2 9l+3 





l = k + 2 

l^k-1 
l = k-2 
other cases. 



(13) 



For irrational Q, the kinetic term Tk has a pseudo-random behavior in k and 
we find that its probability distribution is narrowly peaked as required for 
localization. The hopping term now couples each component of the eigenstate 
with its two closest neighbors. An equation similar to eq. (11) holds for the 
other chirality component. We have shown, therefore, that the modified ran- 
dom walk can be mapped into a localized one-dimensional Anderson model. 



5 Conclusions 



A modified version of the discrete time quantum walk on the line has been 
considered. After an initial diffusive stage, dynamical localization takes place 
for most values of the parameter Q. However, for rational values of the pa- 
rameter the dynamics of the system corresponds to a resonant behavior. In 
particular the usual quantum walk is recovered when Q is an integer. The 
modified walk docs not show a quadratic increase of its variance for all values 
of the parameter; rather, this happens only for exceptional (rational) values, 
for which the quantum spreading is faster than the classical diffusion spread. 

The modified quantum walk has been mapped into a one-dimensional Ander- 
son model, as was previously done in the case of the kicked rotor [24,25]. The 
similarity between the modified quantum walk and the kicked rotor has thus 
been established both numerically and analytically. Both models are charac- 
terized by two parameters, a strength parameter and a scale parameter. In 
the case of the modified quantum walk, the scale parameter is = . 
The strength parameter can be associated to the parameter determining the 
generalized rotation of chirality, as discussed in [15]. 

In the generic case, when Q is irrational, a new kind of quantum random walk 

on the line, with exponentially localized wave functions, is obtained. This 
feature opens interesting possibilities for quantum information processing. In 
particular the question arises, may DL be used for controlling the spread of a 
wavefunction? Furthermore, since DL states are robust against moderate de- 
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coherence [9,23], can they be used in preserving quantum information against 
environment induced noise? It would be interesting to experimentally observe 
the effects of resonances and DL when one of the existing experimental pro- 
posals for the quantum walk [2,3] is realized. We note that these effects might 
also be observed in the context of classical wave simulations for the quantum 
walk [5,6]. 

We acknowledge support from PEDECIBA and PDT S/C/OP/29/84. R.D. 
acknowledges partial financial support from the Brazilian National Research 
Council (CNPq) and FAPERJ (Brazil). A.R. acknowledges useful comments 
made by V. Micenmacher. 
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